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A quantum kinetic theory of the spin transfer between carriers and Mn atoms in a Mn doped 
diluted magnetic semiconductor is presented. It turns out that the typical memory time associated 
with these processes is orders of magnitude shorter than the time scale of the spin transfer. Nev- 
ertheless, Markovian rate equations, which are obtained by neglecting the memory, work well only 
for bulk systems. For quantum wells and wires the quantum kinetic results qualitatively deviate 
from the Markovian limit under certain conditions. Instead of a monotonic decay of an initially pre- 
pared excess electron spin, an overshoot or even coherent oscillations are found. It is demonstrated 
that these features are caused by energetic redistributions of the carriers due to the energy-time 
uncertainty. 

PACS numbers: 75.78.Jp, 75.50.Pp, , 78.47.J-, 75.30.Hx 



I. INTRODUCTION 

Ultrafast dynamics of magnetic semiconductors are 
a vastly growing field. Most studied are manganese 
doped diluted III-V or II- VI semicondutors and their 
nanostructures»ii^ Driven by the vision of spintronic ap- 
plications, much of the research is devoted towards the 
control of magnetic properties on ultrashort time scales. 
Here, optical methods are particularly attractive. Exper- 
iments as well as theoretical studies have revealed that 
optical manipulation of the magnetization is possible in a 
coherent non-thermal regime^^— which opens many new 
perspectives compared with schemes that rely on thermal 
effects onlyJ^^— While most of the coherent magnetiza- 
tion dynamics studied so far in extended semiconductors 
are dealing with the observation and the control of co- 
herent spin precessiou)^^— the transfer of spin between 
the Mn and the electron or hole subsystems is usually 
considered to be an incoherent assimilation process that 
can adequately be described by a Markovian ratej^^rJ^ 
In contrast, for single quantum dots with a single em- 
bedded Mn atom the spin transfer is a coherent process 
where spin is exchanged between discrete levels by Rabi- 
type rotations^^r^ It will be shown in this paper that 
also in quantum wells and wires a coherent spin transfer 
between carriers and Mn atoms is possible where spins 
are exchanged back and forth between these subsystems. 
Our analysis reveals a pivotal role of the energy-time un- 
certainty for enabling the coherent spin-exchange in ex- 
tended systems. 

In this paper, we concentrate on the spin transfer be- 
tween Mn atoms and electrons in the II-Mn-VI semi- 
conductor ZnMnSe and disregard other mechanisms that 
may change the electronic spin. This spin transfer pro- 
cess, of course, in general competes with many other pro- 
cesses that may affect the electronic spin dynamics such 
as the Dyakonov-Perel, Elliott- Yafet or Bir-Aronov-Pikus 



mechanismii which may lead to a spin relaxation. In 
Ref. [l^ it was shown that in the conduction band of II- 
Mn-VI semiconductors the spin transfer can become the 
dominant intrinsic influence on the electron spin which 
justifies our focus on this mechanism. 

We shall demonstrate that the spin transfer may pro- 
ceed qualitatively different from what is predicted by a 
Markovian-rate equation even though the memory asso- 
ciated with this process is short. While a short mem- 
ory time is a necessary precondition for the validity of 
Markovian-rate theories we still find deviations from the 
rate behavior that can be traced back to the fact that 
the exchange interaction which is responsible for the spin 
transfer between electrons and Mn atoms simultaneously 
leads to a redistribution of electron energies that is lost 
in the Markov limit. We have performed simulations 
of the pertinent dynamics for bulk, well, and wire sam- 
ples and find the most pronounced deviations from the 
Markov limit for quasi one-dimensional wire systems. For 
wells the deviations are weaker but clearly visible while 
bulk sample are well described by the Markovian theory. 
We are using for our calculations a quantum kinetic the- 
ory of the correlated spin dynamics in diluted magnetic 
semiconductors (DMS) that we have recently developed 
which treats the exchange interaction between Mn atoms 
and carriers beyond the mean-field and virtual crystal 
approximations.— While other approaches that explicitly 
account for correlations have been derived in the frame- 
work of the Greens-function formalism^ our approach 
is based on a microscopic density matrix theory. 



II. 



MODEL 



In this paper, we apply the theory developed in Ref. 1 13 
to the description of the spin transfer dynamics in the 
conduction band of ZnMnSe. We consider a single band 
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model a Hamiltonian H consisting of two parts: 

H = Ha + H,d, (la) 

the electronic band structure Hq of the host semiconduc- 
tor and the exchange interaction Hsd between s-like con- 
duction band electrons and c?-like localized Mn orbitals. 

The host semiconductor may be a bulk material, a 
quantum well or a quantum wire which, when concentrat- 
ing on the lowest subband, represent systems of effective 
dimension D = 3,2, 1, respectively. In all cases, we ac- 
count for two spin-degenerate conduction bands with a 
well defined spin quantum number of s = ^. In second 
quantization, Hq then reads: 

i^o =^^^k4kC,k, (lb) 

o-k 

where a € {J,, t} is the spin quantum number, k denotes 
the ZJ-dimensional wave- vector and the corresponding 
carrier energy. 

The Mn atoms have half-filled d-shells with vanish- 
ing angular momenta and can be described as localized 
S = I spins in good approximation. The exchange inter- 
action in contact form between these spins and the con- 
duction band electrons is given by the Kondo-like Hamil- 
tonian: 

H- , - V S ,-S ,p»(k'-k)Ri t ^,,6/ 

^^sd — Y / y '-'mi' »crcr'C ''crk f k' rm' i \^^} 

Inn' 

where Jsd is the exchange constant, V the sample Vol- 
ume, S„„' the vector of Mn spin matrices with n S 
{— |, . . . , |}. Sccr' denotes the vector of spin Pauli ma- 
trices and R/ stands for the position vector of the I-th 
Mn atom. The operator 

Pi„, :=|/,n)(/,n'| (2) 

represents the Mn spin degrees of freedom, where |/,n) 
denotes the rt-th eigenstate of the z-component of the 
I-th Mn spin. In real samples the Mn dopants are ran- 
domly placed all over the sample. The Hamiltonian Hsd 
is formulated for a given specific configuration. 

Starting from the Hamiltonian in Eq. ([T]), we have de- 
rived a closed set of quantum kinetic equations of motion 
[Eqs. (16)-(17) in Ref. [l^ which describe correlated Mn 
and electron dynamics in an on average spatially homo- 
geneous systemiiS The averaging involves a quantum me- 
chanical average as well as an average over the random 
Mn positions which are assumed to be homogeneously 
distributed in space. 

In order to focus on the spin transfer dynamics we 
consider a situation where initially all electron spins arc 



aligned along the z axis and the Mn spin are equally dis- 
tributed over the six spin states leading to a zero Mn 
magnetization. With these initial conditions no preces- 
sion will occur as the electron spin vector will stay at 
all times parallel to the z-axis and to the Mn spin vec- 
tor that will acquire finite values in the course of time 
due to the spin transfer. As a consequence, it turns out 
that only the following subset of the dynamic variables 
of Eqs. (16)-(17) in Ref. has non-zero values: 

M: (PL), (3a) 
Ck := (clkCak), (3b) 
C^tr-^Vdicl^^w'^^^-"^^^^'), (3c) 

K::t, V^<5(cI,k,c..k.P„',„,e^('^-''^)^0, (3cl) 

where (...) stands for both the quantum mechan- 
ical average with respect to the statistical op- 
erator and the disorder average with respect to 
the random spatial distribution of the Mn atoms. 
The correlation functions 6{cl^^^Ca-]f^2^^^^^~^^^^') and 

'5(cIikiC(T2k2-Prfi«2e'^'^^~'''^^'> occurring in Eqs. §^ and 
pd[) are related to the corresponding expectation values 

/„t „ „i(k2-ki)R7\ „ J / t pi j:(k2-ki)R7\ 

V'^crki'^'Tkae / ana \C^jki'^'T2k2-f^nin2'^ / 

by subtracting from the latter all possible factorizations 
into factors involving expectation values of fewer opera- 
tors. The lengthy explicit definitions of these correlation 
functions can be found in Eq. (13) of Ref. IT9I . 

Af " represents the average occupation of the n-th Mn 
eigenstate and C^j^ the occupation of the Bloch state 
with wave vector k and spin cr. 

C^k? describes a corre- 
lation between electrons and Mn positions R/ (averaged 
over the random positions of the Mn atoms) which is 
due to the scattering of the electrons with the localized 
Mn atoms. This correlation function would also arise if 
the interaction between electrons and Mn spins would be 
given by a spin independent localized potential at ran- 
dom Mn positions instead of the exchange interaction 
Eq. ([Tc| . K^^^^^'^ describes correlations between elec- 
trons, Mn spins and their respective positions. In this 
function also spin correlations between the Mn and the 
electron subsystems are represented that arise because 
the exchange interaction in addition to providing a local- 
ized scattering channel acts on the spin degree of free- 
dom. 

Using this reduced set of variables the pertinent equa- 
tions of motion can be derived for bulk, well and wire 
systems along the lines detailed in Ref. [l^ resulting for 
the D dimensional case in: 
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(4b) 



/^cr 1 y^cr 1 I 

^<7k2 ^cTiki y(27r)-^ 



BZ 



BZ 



(4c) 
(4d) 



where the source Q ^^2^2^^ on the right hand side of 

Eq. (j4dp is given exphcitly in the appendix. The inte- 
grals are to be taken over the first Brillouin Zone {BZ). 
L denotes the length of the system in the unconfined di- 
rections, V ~ for bulk systems, V = L^d for quantum 
wells of the thickness d and V = LA for quantum wires 
with a cross sectional area A. For the numerical imple- 
mentation of the above equations it is advantageous to 
note that in the case considered here, K is only non-zero 
if the indices a and n fulfill either 

ci = o'2, fii = 71,2 or, (5a) 

CTl = (72 ± 1, ni = 7l2 T 1, rii + (Ti = 712 + (72. (5b) 



III. GENERAL PROPERTIES OF THE 
EQUATIONS OF MOTION 

A. Relative importance of different correlations 

In this subsection we shall demonstrate that for the 
spin transfer processes to be discussed in this paper only 
a much smaller subset of terms essentially contributes. 
The resulting reduction in complexity is not only advan- 
tageous for the numerics, it will also enable us to analyze 
more conclusively the pertinent features of the resulting 
dynamics. 

The key observation for identifying the most important 
contributions is that for an initially uncorrelated system 
in which the average Mn spin is zero and the electrons are 
completely spin polarized the onset of the spin transfer is 
mediated exclusively by the first source term in Eq. (j4d[) 
for the correlation function K 
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(6) 



as has been shown in Ref. jl9|. This implies that this term 
plays a pivotal role for the spin transfer dynamics that is 
the target of the present paper. The resulting K directly 
drives via Eqs. ([^a]) and (|4b[) the main observables of in- 
terest M and C. The rather involved additional sources 
Qf^ (cf. the appendix) in Eq. (j4dp as well as the dis- 
order related correlations C build up subsequently in a 
second step and are of higher order in the coupling con- 
stant Jsd^ Based on this observation it is suggestive to 
assume that Qj^ and C might be of less importance also 
when the dynamics is followed over a longer time scale, 
although the correlation functions and C appear on 
the same level of the correlation expansion as the remain- 
ing contributions. In order to test this assumption, we 
have performed numerical simulations where we compare 
results of the full set of equations Eqs. Q with calcula- 
tions where and C have been set to zero. 

To be specific, we have numerically solved the initial 
value problem of initially spin polarized electrons with a 
Gaussian distribution in energy space according to 
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(7a) 



and Mn spins where initially all possible z components 
have equal probabilities, i.e. 



= g- (7b) 
Finally, all correlations are initially set to zero 

C:t\to = KZtk - 0. (7c) 



These initial conditions are meant to mimic the situation 
immediately after a fast optical excitation of an initially 
unmagnetized DMS. 
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FIG. 1. (color online) Time evolution of the total electron spin in Zno.93Mno.07Se for (a) a bulk system; (b) a quantum well; 
(c) a quantum wire. The solid red lines have been calculated with the full quantum kinetic Eqs. Q, the blue lines with bullets 
have been obtained by neglecting and C in Eqs. Q , and the green lines with squares represent the results of the Markovian 
rate Eqs. I|14p . The pink curve with triangles in (b) has been obtained from Eqs. (|19|l . where the energetic redistribution of 
electrons is neglected. 



We have performed simulations for a wide range of pa- 
rameters typical for DMS. In particular, we have varied 
the exchange coupling Jsd in the range 5-100 meVnm'^ 
and the Mn concentration x in the range of 1-10 % and 
used an effective electron mass of me = 0.21mo with 
mo being the free electron mass, which is typical for 
Zni_j;Mn3;Sc. The parameters A and Eq of the initial 
Gaussian distribution in Eq. (|7ap have been taken to be 
A = 0.4 meV and Eq = meV, respectively. The above 
initial value problem has been solved for bulk, a well with 
width d = 4 nm and a wire with a cross-section of area 
A = 16 nm^. It turns out that the results are qual- 
itatively similar in the whole parameter range studied. 
Typical results that correspond to Zno.93Mno.07Se where 
Jsd = 12 meVnm'^ are shown in Fig. [U where the total 
electron spin 

Se=^S,,Ck (8) 
kcr 

is plotted as a function of time for (a) a bulk semicon- 
ductor, (b) a quantum well and (c) a quantum wire. The 
solid red lines represent results of the full set of equations 
while the blue thin lines with bullets are obtained by ne- 
glecting Qj^ and C. The green lines with squares are 



obtained from Markovian rate equations and will be dis- 
cussed later. At this point we only want to note that the 
results of the former two calculations quantitatively are 
almost indistinguishable for all conditions studied. For 
a detailed discussion of the physical implications of the 
curves in Fig. [1] it turns out to be useful to exploit the 
fact that Qj^ and C can be safely neglected in our case 
for a reformulation of the remaining equations. 

B. Integral representation and memory kernel 

When Qj^ and C are discarded the right hand side of 
the equation of motion (j4d[) for K is independent of K 
and can easily be integrated resulting in: 

X (^^^^k. (^')M„",^ it') Cll^^ {t')A'Q it')) dt'. (9) 

Here we have again assumed that A'^^^"^^^ (to) — 0. In- 
serting Eq. ([9]) in the equations of motion (|4a| and (j4bp 
for C and M formally eliminates K from these equations 
leading to a closed set of equations of motion involving 
only C and M as dynamical variables: 
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BZto 



(10b) 



9t '"I'^i hWi2TT)i^ 

BZto 



where the k-dependent memory function (r) is de- 
fined as: 

G^(t) := I cos (i(i?k' - i?k)r)rf^k'. (11) 

BZ 

If one is interested in the total electron spin S*^ in 
the system Eq. (jlObp has to be summed over k. After 
this summation the same memory function (r) as in 
Eq. (jlOap appears also in the equation for S*^ demon- 
strating that G^ (r) determines the memory of the spin 
transfer process. Formulating the equations of motion us- 
ing a memory kernel has the advantage that we can now 
analyze the corresponding memory depth. In the case of 
a quantum well this can be made even more explicit as 
for parabolic bands (i.e. £'k = ) ^^'^ approximating 
the Brillouin zone as a sphere with radius ksz — Iksz] 
the integral in Eq. pT|) can be evaluated analytically re- 
sulting in: 



GUt) 



BZ 



2^ 



hT). (12) 



Fig. [2] shows the memory kernel G1{t) as a function of r 
using the ZnMnSe effective mass me = 0.21too for three 
representative values of k [k = {k — ksz/V^) cor- 
respond to the k values where the first zero of G^(t) 
comes at the earliest (latest) time while k = OAksz is 
an intermediate value] . In all cases G^ (r) decays and os- 
cillates rapidly on a sub-femtosecond timescale and de- 
pends only weakly on k. These general statements also 
hold for bulk and quantum wire systems, where we have 
calculated G^ (r) numerically. 

As can be seen in FigUl the timescale of the exchange 
induced spin dynamics in typical DMS is of the order of 
several picoseconds, i.e. over 10^ times larger than the 
memory depth of the system. One could, thus, suspect 
that it is justified to treat the spin transfer dynamics in 
the Markovian limit. 
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FIG. 2. (color online) Spin memory function G|(r) of a Zn- 
MnSe quantum well for different values of k. 



the results with our quantum kinetic theory represented 
by Eqs. (j4|) that here give essentially the same results as 
the non-Markovian Eqs. (cf. Fig.p. 

The Markov limit corresponds to neglecting the retar- 
dations in Eq. ([9]) of the functions M and G. Formally, 
this is done by evaluating M and C at time t instead of 
t' which allows to take these functions out of the integral 
over the past. Finally, the limit to ~> ~oo is taken which 
eliminates the initial time. The remaining integral yields 
(5-distributions and principal value integrals. The latter 
do however drop out when inserted in the equations of 
motion for M and C. The principal parts also vanish 
in the limit i — > oo which is the limit usually discussed 
when deriving Fermi's golden rule. For a parabolic band- 
structure this procedure yields in Eq. ([9]): 



j>^cr2n2k2 



X inS 



2m, 



(13) 



C. Markovian limit 

In this subsection we shall derive the Markovian limit 
of the spin transfer described by our model and compare 



It should be noted that going over to the Markov limit in- 
troduces the energy conserving delta function in Eq. (|13p 
which for parabolic bands makes K diagonal with respect 
to the absolute values of ki and k2. 

With K approximated by Eq. the equations of 
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motion (|10ap and (jlObp for M and C then become: 
d 



dt "1 



J2 m,r!z3L2i5 



S3F2(27r)2-D- 



(14a) 



BZ 



nn' (7 



fT 1 A; 1 



D-2 



(14b) 



where = 47r, = Stt and fii — 2. In the above equa- 
tion, we have assumed the bandstructure to be parabohc 
and averaged the electronic variables with respect to the 
angle of their k- vectors: 



ak — 



n, 



For such a bandstructure, this averaging procedure di- 
rectly leads to a closed set of equations of motion without 
requiring any further approximations. Since the num- 
ber of Mn atoms in an DMS is usually several orders of 
magnitude larger than the number of optically generated 
spin polarized electrons, the Mn variables are almost 
constanti^ In the paramagnetic case with zero magnetic 
field, M"(t) is approximately i. Then, it follows from 
Eq. (jl4bp that the spin of an electron with momentum 
hki is given by s^^ = s|^C|j^^^ sl^Cf^ 
exponentially: 



^^^^ simply decays 



with the rate 



7fci 



D- 



sd' 



-rtMnki 



D-2 



(15) 



(16) 



This rate, which we have obtained as a limiting case of 
our full quantum kinetic equations ([4]), is already known 
since the 1970s where it has been derived from a golden 
rule analysis for the bulk casci^ For a quantum well the 
rate in Eq. p6)) has been deduced in, e.g., Ref. [13. 

Now, the quantum kinetic Eqs. ([4]) can be used to check 
the validity of the rate Eqs. ([H]) for the description of 
spin transfer dynamics in different DMS systems. Fig. [T] 
shows the electron spin dynamics in Zni_2;Mn3;Se sys- 
tems of different dimensionality and compares quantum 
kinetic results (thick solid lines) with the Markov limit 
represented by the rate Eqs. (|14p (green dashes lines with 
squares). As can be seen: (a) the rate equations are 
in excellent agreement with quantum kinetic equations 
for bulk semiconductors, (b) for a quantum well the full 
equations predict a clear overshoot of the electron spin 
which is absent in the Markovian limit and (c) accord- 
ing to the full equations, the electron spin in a quantum 
wire can even show oscillations that are also not expected 
from rate equations. 



Before analyzing in detail the origin of these non- 
Markovian features, we shall shortly examine the impact 
of material parameters. If we again approximate M as 
constant and rescale the time t with mo /me, it can be 
seen from Eq. (jlObp that material parameters enter the 
electron spin dynamics only as the prefactor 



(17) 



of the source terms on the right hand side. Here, 
we illustrate the influence of material parameters ex- 
emplarily by simulating the initial value problem from 
Eqs. ([7]) for quantum wells with different prefactors F = 
meJg^nMnd~^ but fixed mass ratio mo /me- Displayed in 
Fig. [3] is the resulting total electron spin as a function of 
the unsealed time t. It can be clearly seen that the am- 
plitude of the overshoot non-monotonically depends on 
F, while the time it takes until the maximal overshoot is 
reached decreases with F. For small values of < O.Si^o, 
practically no visible overshoot takes place and the re- 
sults are in good agreement with the rate equations (jl4p , 
for F « I.8F0 the overshoot amplitude reaches its max- 
imum before it decreases with increasing F. A similar 
trend is found for quantum wires, where the strength 
of the oscillations depends non-monotonically on F (not 
shown) . 




t (ps) 

FIG. 3. (color online) Time evolution of the total electron spin 
for different values of the rescaled electron Mn coupling F 
[cf. Eq. (|17[l ] at fixed mass ratio mo/me. Fq corresponds to 
the parameters used in Fig. [T](b). 

The spin overshoot and oscillations are clear signatures 
of coherent dynamics which cannot be described on the 
level of rate equations. It should be recalled that the 
oscillations are not related to a precession of the spins 
which is not possible for the configuration considered 
here as the Mn and electron spins are aligned parallel. 
Instead they represent a coherent exchange of spin be- 
tween the electronic and the Mn subsystems. It is also 
worth noting that the rate equations fail in the two- and 
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one-dimensional case even despite the femtosecond spin 
memory depth revealed by Eq. (|lip and seen in Fig. [2] 
for D = 2. 



IV. ORIGIN OF THE SPIN OVERSHOOT 

As a short memory time is commonly believed to guar- 
antee the validity of the Markov approximation the ques- 
tion arises, why the rate equations fail to describe the 
dynamics in the two- and one-dimensional systems al- 
though this condition is fulfilled. It turns out that the 
corresponding analysis yields qualitatively the same an- 
swers for wells and wires. Here we shall discuss explicitly 
only the case of quantum wells as these systems are more 
widespread while Mn doped quantum wires are still a 
very novel field of researchi^ 

The key to understand the origin of the observed non- 
Markovian behavior is to analyze not only the total elec- 
tron spin Se, which is obtained as a summation over k 
space [cf. Eq. ©J but to follow the time evolution of the 
electron occupation and the electron spin over the kinetic 
energies £^k of the electrons. 

Fig. |4] shows the electron distribution as a function of 
iJk and the time t for the parameters used for Fig. [T](b). 
Initially, the electrons are Gaussian distributed. They 
are subsequently scattered towards higher energies, fol- 
lowed by a reflux. These tails occur repeatedly before a 
quasi static electron distribution is reached for t > 25 ps. 
Eventually 19 % of the electrons in the system have been 
redistributed within the Brillouin zone. 

This energetic redistribution is due to the scattering of 
electrons with the spatially localized potentials of the Mn 
atoms. As discussed in Ref. [l^ even in the limit of a spa- 
tially homogeneous Mn distribution which result in an on 
average spatially homogeneous system k-space scattering 
takes place as a result of the localization of the Mn scat- 
tering centers. This scattering has been shown to be ac- 
companied by a redistribution over electronic energiesi^ 
which in a quantum kinetic description is allowed due to 
the energy-time uncertainty. We note that, according to 
the rate equation (|14bp , no such electron redistribution is 
predicted since in the Markovian limit the kinetic energy 
is conserved in the scattering process [cf. Eq. P^ ]. 

The electron spin distribution in energy space that cor- 
responds to the electron distribution in Fig. 2] is depicted 
in Fig.[5j It does not simply decay exponentially as would 
be expected in the Markovian limit according to Eq. ([T5]) . 
Instead, the spin distribution forms several tails in the 
course of time which correspond one-to-one with similar 
tails of the electron occupations in Fig. |4l Each tail can 
be associated either with a positive or negative electron 
spin. The sign of the spin of consecutive tails alternates, 
with the spin of the first tail being opposite to the positive 
initial spin. This can be explained by the nature of the in- 
teraction Hsd of the electrons with the scattering centers, 
i.e. with the localized Mn spins. According to Eq. (ITcl) . 
the scattering of the itinerant electrons in k space can 




FIG. 4. (color online) Time evolution of the electron distri- 
bution in energy space. 
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FIG. 5. (color online) Time evolution of the electron spin 
distribution in energy space. The electron spin is given in 
multiples of |. 



take place with or without a flip of their spins due to the 
spin exchange with the Mn atom. Analyzing the struc- 
ture of the involved spin matrices reveals that processes 
with spin flip occur more often than those without spin 
flip. Hence, starting initially with a positive spin the first 
tail resulting from the scattering is predominantly nega- 
tive. The second tail contains mainly electrons that have 
been rescattered from the first tail to higher energies and 
consequently their spin is positive. 

In order to analyze these redistributions in more detail 
we have plotted in Fig. [5] (a) the electron and in Fig. [S] 
(b) the spin distribution for different times t, i.e. cross 
sections of Figs. |4] and [SJ respectively. While the initial 
electrons are Gaussian distributed [cf. Eq. ([7a|) ]. the fi- 
nal electron distribution is a Lorentzian curve in good 
approximation. The local minima and maxima in the 
electron occupation, which can be clearly identified in 
the inset of Fig. |6] (a), correspond to the tail structure 
in Fig. m Although the total electron spin is basically 
gone after 20 ps [cf. Fig. [T] (b)], the spin distribution is 
not zero yet for all energies but rather quickly oscillates 
with respect to the energy [cf. inset in Fig. [S^b)]. With 
increasing t, these oscillations become fast and their am- 
plitude slowly vanishes. Thus, eventually the electron 
spin distribution approaches zero which implies that the 
initial electronic spin is completely transferred to the Mn 
subsystem. 
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FIG. 6. (color online) Distribution of (a) the electron occupation and (b) the electron spin as a function of the kinetic energy 
for different times t. The insets show magnified sections of the main figures. 



This behavior should be contrasted with spin assim- 
ilation processes described by rate equations where the 
electron energy is conserved. In such a process the major- 
ity (minority) spin at a given energy decreases (increases) 
until any excess spin is gone. Thus, at each energy the 
total spin cannot switch its sign. In contrast, quantum 
kinetics allows that the majority spin at a given energy is 
flipped and transferred to a different energy. Thus in cer- 
tain regions of energy majority and minority spins switch 
roles repeatedly and it depends on details of the distri- 
bution whether or not the total electron spin switches 
sign. Indeed, the overshoot in quantum films and the 
oscillations in quantum wires demonstrate that such a 
reversal of the sign of the total electron spin may actu- 
ally take place [cf. Fig. [Tj. These observations suggest 
that the most severe approximation implied by replacing 
the quantum kinetic memory function by a Markovian 
rate is not the neglect of a finite memory time but the 
suppression of the redistribution in energy space that in 
the quantum kinetic theory accompanies the electron-Mn 
scattering. In order to substantiate this assumption we 
have developed a level of description where the memory 
depth is kept finite but the energy redistribution is arti- 
ficially suppressed. To this end we make the ansatz 



(18) 



for the electron variables where we force the electron dis- 
tribution /(k) to be time independent which implies that 
the energy distribution is kept fixed. Without loss of gen- 
erality /(k) can be assumed to be normalized according 
to 



1= y'/(k)d^k. 

BZ 



Eq. (jlObp over ki wc obtain 
_M"i - '''' V (S 



c^{t')M:{t') - c::{t')M:i{t'))G''{t - t')df 



to 



(19a) 



d_~ 
dt 



" h^V(2TT)D J 2^ ^^^') 

to 

X (c^m: - c:im:;) c^it - t')dt'. (1%) 

Here we have introduced the k-integrated memory func- 
tion 



G^(r) 



G^(r)/(k)d^k, 



(20) 



BZ 



Inserting this ansatz in Eqs. (jlO[) and integrating 



where (r) was defined in Eq. pT|) . Solving these equa- 
tions for the example of a quantum well with the initial 
conditions Eqs. ([7|) [cf. pink line with triangles in Fig. [T] 
(b)], essentially reproduces the results, we have previ- 
ously obtained from the corresponding rate equations, 
even though we now account for the finite memory depth 
and only neglect the redistribution of electron energies. 
In particular, the spin overshoot can only be described 
with the full quantum kinetic theory. It thus has to be 
concluded that these redistributions are indeed the ac- 
tual source of the reversal of the sign of the electron spin 
observed in the quantum kinetic simulations and that 
the neglect of the finite memory time is the less severe 
approximation when deriving the Markovian rate equa- 
tions. 
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t (ps) kinetic energy (meV) kinetic energy (meV) 

FIG. 7. (color online) (a) Time evolution of the total electron spin for different central energies Eq of the initial Gaussian 
distribution; (b) distribution of the electron occupation as a function of the kinetic energy for the central energy Eo = 3 meV 
and different times t; (c) the corresponding distribution of the electron spin 



The same analysis for bulk and quantum wire systems 
reveals that energy redistribution of electrons and corre- 
sponding spin oscillations in the Brillouin zone are also 
present, but they are far less pronounced in the 3D case. 
Hence, the rate equations agree with the full quantum 
kinetic equations [cf. Fig. [T] (a)] for bulk systems and 
predict a monotonic spin assimilation in this case. In 
quantum wire systems, on the other hand, these effects 
are much stronger, leading to an oscillating total elec- 
tron spin [cf. Fig. [T] (c)]. It is interesting to compare 
our present results also with the Rabi-type spin exchange 
observed in single quantum dots doped with a single Mn 
atomii"— In a Rabi-type spin transfer essentially two dis- 
crete states are involved that are coupled by the exchange 
interaction. The dynamics of this two-level system can 
be reduced to the dynamics of the occupations by for- 
mally solving the equation of motion for the coherence 
between these states. Inserting the result into the equa- 
tions of motion for the occupations yields a memory func- 
tion representing the coherence. Noting, that in our case 
where the Mn concentration is much larger than the itin- 
erant electron concentration and therefore changes of 
are negligible in determining the back action of the Mn 
atoms on the carriers, we find according to Eq. (jlObp that 
also in our case the change of the electronic occupations 
at time t is determined essentially only by their values 
at earlier times. However, in the two-level case there 
are only two occupations and these are not independent 
variables, because of the charge conservation. Thus, in 
this case a given occupation is coupled only to its own 
values at earlier times. Therefore, the non-monotonic 
time evolution of the electronic spin as manifested in the 
corresponding Rabi-flops can only take place when the 
memory provided by the coherence is long, as otherwise 
we would necessarily reach the limit of Markovian rates. 
An energetic redistribution of carriers, as required for co- 
herent spin transfer in quantum wells and wires, is not 
possible because there are no further final states available 
in a two-level system. We therefore conclude that the 



mechanisms responsible for coherent spin exchange differ 
in extended semiconductors qualitatively from those in 
quantum dots. 

We note in passing that we have performed a series of 
simulations for different shapes and widths of the initial 
distributions [not shown] which lead to the conclusion 
that, as a general rule, the redistribution of electronic 
energies and the resulting effect on the spin dynamics 
is more significant the sharper occupied and unoccupied 
regions are initially separated in the Brillouin zone. For 
example, relatively broad initial distributions which re- 
semble a Lorentzian curve do not tend to show a spin 
overshoot, while the overshoot is especially pronounced 
for narrow box like distributions with steep edges. 

Apart from the shape and width another important 
parameter of the initial distribution is its energetic posi- 
tion determined by the central energy Eq of the Gaussian 
electron distribution given in Eq. (ffa)) . Fig. [7] (a) shows 
the total electron spin as a function of time for different 
values of Eq. The curve with Eq = nieV has already 
been shown in Fig. [T] (b) and is repeated here for better 
comparison. For Eq — 0.4 mcV, which equals the stan- 
dard deviation A of the Gaussian distribution, the spin 
overshoot is less distinct. For £'o ^ A, the time evolution 
of the total electron spin no longer depends on Eq and 
the overshoot disappears. The electron and spin distribu- 
tion for Eq = 3 meV are depicted in Figs. [7] (b) and (c), 
respectively. It is seen that energetic redistributions take 
place also in this case and the spin distribution oscillates 
as a function of time and energy. Compared with the 
corresponding distributions for _Bo = meV [cf. Figs. [S] 
(a) and (b)] these distributions have a higher symmetry 
and the excess spin initially prepared at the central en- 
ergy Eq decreases faster as now redistributions to lower 
and higher energies are possible. It turns out that these 
quantitative differences altogether have the effect that 
here after summing over the individual spins the oscilla- 
tions average out and the total spin simply decays mono- 
tonically. A similar behavior is found in quantum well 
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systems, where the osciUations of the total electron spin 
are most pronounced if the initial electron distribution is 
in the vicinity of the band edge (not shown) . This obser- 
vation suggests that the different density of states at the 
band edge in bulk systems, quantum wells, and quantum 
wires is is an important factor for the appearance and 
strength of the coherent phenomena in the spin transfer. 

V. CONCLUSION 

We have analyzed the spin transfer from an initially 
prepared electronic excess spin towards the spin of Mn 
atoms in a diluted Mn doped magnetic semiconductor 
within a quantum kinetic theory. We have demonstrated 
that these spin transfer processes are dominated by only 
a few relevant correlations which form a small subset of 
the full set of terms that contribute to the quantum ki- 
netic equations on this level of the correlation expansion. 
By concentrating on these terms only, it is possible to 
eliminate the pertinent correlations in favor of a memory 
function. Although the typical memory times are orders 
of magnitude shorter than commonly encountered spin 
transfer times, noticeable deviations from Markovian dy- 



namics may occur for quantum well and wire systems 
while for bulk systems we find a good agreement be- 
tween quantum kinetic calculations and the Markovian 
limit. For wells and wires the spin transfer involves one 
or more changes of the sign of the total electron spin 
which indicates a regime of coherent exchange of spin be- 
tween the electronic and the Mn subsystems. In contrast 
to Rabi-type spin-exchange in Mn doped quantum dots, 
these coherent dynamical features are here not related 
to a long memory time. Instead, it is a redistribution 
of electronic energies due to the energy-time uncertainty 
which enables a spin-exchange where the total electron 
spin exhibits a non-monotonic time evolution. Marko- 
vian rate equations fail to describe this type of dynamics 
because the energy-time uncertainty is suppressed in this 
limit and thus the necessary energetic redistributions do 
not take place. 

Appendix A: Source terms for correlations 

In this appendix we give the explicit form of the source 
term on the right hand side of the equations of mo- 
tion (|4d[) for the correlations K. It reads: 
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